FROM A KINETIC EQUATION TO A DIFFUSION 
UNDER AN ANOMALOUS SCALING 

GIADA BASILE 

Abstract. A linear Boltzmann equation is interpreted as the 
f>o ", forward equation for the probability density of a Markov pro- 

cess (K(t),i(t),Y(t)) on (T 2 x {1,2} x R 2 ), where T 2 is the two- 
^^ , dimensional torus. Here (K(t),i{t)) is an autonomous reversible 

jump process, with waiting times between two jumps with finite ex- 
£^ . pectation value but infinite variance. Y(t) is an additive functional 

of K, defined as L v(K(s))ds, where \v\ ~ 1 for small k . We prove 
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that the rescaled process (N In N)~ 1 / 2 Y(Nt) converges in distri- 
bution to a two-dimensional Brownian motion. As a consequence, 
the appropriately rescaled solution of the Boltzmann equation con- 
verges to a diffusion equation. 



1. Introduction 



One of the most interesting aspects of the problem of energy trans- 
port in a solid is the anomalously large thermal conductivity observed 
in low dimensional materials (see [20] , [7] for a general review; see also 
[TT] for experimental data for graphene materials). So far very few re- 
sults are obtained by a rigorous analysis of microscopic dynamics, and 
even crucial points, such as the exponent of the divergence of thermal 
q ■ conductivity in dimension one, are still debated. 

The theoretical approach proposed by Peierls [25] intended to com- 
pute thermal conductivity in analogy with the kinetic theory of gases, 
conforming to the idea is that at low temperatures the lattice vibra- 
tions, responsible of energy transport, can be described as a gas of inter- 
acting particles (phonons). The time-dependent distribution function 
of phonons solves a Boltzmann type equation, and an explicit expres- 
sion for the thermal conductivity is obtained, which is of the form of 
the kinetic theory n = J dkCkV^Tk. Here Ck is the heat capacity of 
phonons with wave number k, Vk is their velocity and r^ is the aver- 
age time between two collisions. A goal of the kinetic approach is the 
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prediction that the mean free path A^ = v^t^ and thus thermal con- 
ductivity are infinite in dimension one when the phonon momentum is 
conserved. 

Over the last years, several papers are devoted to achieve phononic 
Boltzmann-type equations from microscopic dynamics (see [28J for main 
ideas and tools). In [2], [22], [19] [26] a kinetic limit is performed for 
chains of an-harmonic oscillators, and in [21] a linear Boltzmann equa- 
tion is rigorously derived for the harmonic chain of oscillators with 
random masses. In [5] the authors consider a system of harmonic os- 
cillators in d dimensions, perturbed by a weak conservative stochastic 
noise. The following linear Boltzmann-type equation is deduced for 
the energy density distribution, over the space M. d , of the phonons, 
characterized by a vector valued wave-number k e T d (<i-dimensional 
torus) 

d t u a (t, r,k) + v(k) ■ Vw Q (£, r, k) 
(1) =T^t'E/ dk'R(k,k')[u p (t,r,k')-u a (t,r,k)}, 

a — 1, .., d, d > 2. Equation in dimension one is similar, except for the 
mixing of the components. The kernel R is not negative and symmet- 
ric. Despite the exact expressions of R and v (the velocity), the crucial 
features are that v is finite for small k, i.e. |u| — > 1 as \k\ — > 0, while 
R behaves like \k\ 2 for small k, and like \k'\ 2 for small k'. Naively, it 
means that phonons with small wave numbers travel with finite veloc- 
ity, but they have low probability to be scattered, thus one expects that 
the their mean free paths have a macroscopic length (ballistic trans- 
port). This is in accordance with rigorous results showing that thermal 
conductivity is infinite in dimension one and two ([5], [1]). 

A probabilistic interpretation of (pQ) provides an exact statement of 
that intuition. The equation describes the evolution of the probability 
density of a Markov process (K(t),i{t),Y(t)) on (T d x {1, ..,d} x R d ), 
where (K(t),i(t)) is a reversible jump process and Y(t) is a vector- 
valued additive functional of K, namely Y(t) = J Q ds v(K s ). K and 
i can be interpreted, respectively, as the wave number and the "po- 
larization" of a phonon, while Y(t) denotes its position. In order to 
investigate the property of the process Y(t), one can look at the Markov 
chain {Xj} on T d given by the sequence of states visited by K(t), and 
at the waiting times {r(JQ)}, where t(JQ) is the (random) time that 
the process spends at the i-th visited state. The vector-valued func- 
tion S n = Yli=i T {Xi)v(Xi) gives the value of Y at the time of the n-th 
jump T n = YH=\ T i.^i)i then Y(t) is just the piecewise interpolation of 
S n at the random times T n . 

The behavior of the rate R implies that the stationary distribution 
of the chain is of the form n(dk) ~ \k\ 2 dk for k small, and since the 
average of r(k) goes like \k\~ 2 for fc < 1, the tail distribution of the 
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random variables (r(Xj)v(Xj)} behaves like 

(2) *[\T(X t )v(X t )\>\]~-}-j W>1. 

A 2 

Therefore, in dimension one and two the variables t{Xj)v{Xj) have 
infinite variance with respect to the stationary measure. We remark 
that the variance has the same expression of the thermal conductivity 
obtained in [5]. 

The one dimensional case is discussed in [3] , where the authors prove 
that the rescaled process N~ 2 ^Y(N-) converges in distribution to a 
symmetric Levy process, stable with index 3/2. Convergence of finite 
dimensional marginals has been proven earlier in [16] . Here we consider 
the other critical case d = 2. S n is now a sum of variables with tail dis- 
tribution ~ -h, which means that if they were independent, they would 
be in the domain of attraction of a multivariate normal distribution. 
Looking at the behavior of the variance 
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i T ( X i)MXi)) 2 l {lT(Xi)VaiXt)l <^x } ~lnA, a €{1,2} 



it turns out that the proper scaling contains an extra factor (Inn) 1 / 2 . 
The rescaled process (nlan)~ 1 ' 2 S n t has a central part, given by the 
sum of truncated variables T(Xi)v a (Xi)ln r fx i )v a (x i )\<^/n}i with finite 
variance and an extremal part that goes to zero in probability, due to 
the extra term (Inn) -1 / 2 . This is a standard argument used for sums 
of i.i.d. random variables with tail distribution (J2J), introduced for the 
first time by Kolmogorov and Gnedenko in [T5], that we adapt to the 
case of dependent variables. 

Then we are reduced to the problem of convergence of a sum of 
centered, dependent, bounded random variables to a Wiener process. 
We propose two different approaches. In Section 15.11 we will use an 
abstract theorem due to Durrett and Resnick [8], based on the invari- 
ance principle for martingale difference arrays with bounded variables 
(Freedman, [T3] and (14J). together with a random change of time (see, 
for example, Helland [IB] and Billingsley [6]). The underlying central 
limit theorem for martingale difference arrays can be found in Dvoret- 
zky [9] , [10] (see also [23] , [18] and references therein) . The alternative 
proof, in Section [6j is based on the convergence of the moments to the 
moments of a Brownian motion, under some asymptotic factorization 
conditions, and it uses combinatorial techniques. In this case we will 
only show convergence of the finite dimensional marginals. The multi- 
dimensional generalization is based a Cramer- Wold argument (see for 
example [61, [TJ, [27], [IB]). 

Convergence of {nhin)~ l / 2 S n . to a two-dimensional Wiener process 
is in the Skorokhod Ji-topology. Moreover, since the random times T n 
are sums of positive variables with finite expectation, one can prove, 
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using the arguments in [3], that {nhin)~ l l 2 Y(n-) converges to a two 
dimensional Wiener process in the uniform topology. 

Finally we show that the properly rescaled solution of the linear 
Boltzmann equation in dimension two converges to a diffusion equation. 
The diffusion coefficient is given by an infrared regularization of the 
thermal conductivity obtained in [I], [5], with a proper renormalization 
f lT3|) . Convergence of solutions of linear kinetic equations to a diffusion 
under an anomalous scaling was also proved by Mellet et al [21], using 
an analytical approach. 

We remark that the case d > 3 can be easily treated with the same 
strategy. In particular the rescaled solution of the Boltzmann equation 
converges to a diffusion equation, with a diffusion coefficient given by 
the thermal conductivity obtained in [I], [5]. 

Acknowledgements. I would like to thank Anton Bovier and Nicola 
Kistler for valuable discussions. 



2. The model 
We consider equation ([[]) in dimension two, namely 
d t u a (t, r,k) + v(k) ■ Vu a {t, r, k) 
^ * = 2_\ / dk'R{k, k')[up(t, r, k!) — u a {t, r, k)\, 

Wa = 1, 2, t > 0, x G M 2 , k G T 2 , with a (vector valued) velocity v and 
a scattering kernel R given by: 

(a\ (i\ sin(7rfc a ) cos(7r£; a ) 2 

(4) v a (k) = -^ Vfc G T , Va G {1,2} 

(E^iSin 2 ^)) 

2 

(5) R(k,k') = 16^sin 2 (7rA; a )sin 2 (7rA4), V^'eT 2 . 

a=l 

We denote with (A(t), i(t)) the jump process with values in T 2 x {1, 2}, 
defined by the generator 

(6) Cf(a, k) = V / dk' R(k, k') [/(/5, k') - f(a, k)} , 

with /:{l,2}xT 2 4l continuous on T 2 . The process waits in the 
state (k, i) an exponential random time r with parameter $(fc,i) 



2 „ 2 

(7) $(£;, j) = V(l - ^j) / t/fc' i?(A;, fc') = 8 V sin 2 (7rA; t 



), 



FROM A KINETIC EQUATION TO A DIFFUSION UNDER AN ANOMALOUS SCALING 

then it jumps to another state (J, k') with probability v [i, k;j, dk'] = 
(1-Sij) P(k,dk'), where 

(8) P(k, dk') := c^)- 1 ^, k')dk' = ^^7^*' 

E/3 sln WJ 

Observe that the two processes K(t) and i(t) are independent. Dis- 
regarding the time, the stochastic sequence {X n } n > of states visited 
by K(t) is a Markov chain with value in T 2 , with probability kernel 
P(k,dk'). Since P is regular and strictly positive and defined on a 
compact set, it is ergodic, i.e. there exists a strictly positive probabil- 
ity distribution tt such that P m (k, I) — » n(I) >0,mf oo, for each / 
open interval (see [TIj). Here P m , m > 2, denote the m-th convolution 
integral of P. By direct computation ir(dk) = ^<&(k)dk. 

The process Y(t), with value in M 2 , is an additive functional of K(t) 

(9) Y(t) = F(0) + / rfs v(K s )ds. 

Jo 

In order to investigate its properties, we define two functions of the 
Markov chain {X n } n > , the clock, T n , with values in M + and the posi- 
tion, S n , with values in M 2 

n—l n—1 

T n = Y,zi Hx e )-\ s n = Y,^ <Xi) ®(x e )-\ 

1=0 1=0 

Here {ee}i>o are i.i-d. exponential random variables with parame- 
ter 1, and we take Sq = 0. The clock T n is the time of the n-the 
jump of the process K(t) and it is a sum of positive random vari- 
ables with finite expectation with respect to the invariant measure, i.e. 
E^ei "^(Xx) -1 ] = 1. S n is a two-components vector which gives the 
value of Y(t) at time T n , i.e. S n = Y(T n ). It is a sum of centered 
random vectors whose components show a tail behavior given in (J2J). 
Moreover, the covariance matrix of each of these vectors is diagonal. 
By denoting with T~ l the right-continuous inverse function of T n , i.e. 
T _1 (t) := inf{n : T n > t}, we can represent process Y(t) as follows: 

Y(t) = Sl_ T -i(t)-ij + v(X\r-i(t)-i])(t - T\r-i(t)-l}), 

where ['J denotes the lower integer part. In particular, Y(t) is the 
(vector valued) function defined by linear interpolation between its 
values S n at the random points T n . 

3. Main results. 
For every N > 2, t > 0, we define the rescaled processes 

(10) T N (t) = ±T [m , TZ 1 {t) = ± : T- 1 (Nt), 

( n ) Z ^°) = 7^ S ^ + ( - Nt -\- Nt ^7^W v ( X l N ^)- 
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Observe that Z^ is a two-dimensional continuous vector denned by 
linear interpolation between its values j= S n at the points n/N. 

We assume that the initial distribution /i of the process K t satisfies 
the condition 



(12) / dfj,(k)\k\- 2 < oo, 

which guarantees that E M [eo0(X o ) _1 ] < oo. In particular, it can be a 
delta distribution Sk (dk), with k e T 2 /{0}. 
Let us denote with 



(13) a 2 := lim E„ 



ei^i(Xi N 



Hx 1 



1 !\srni2Sil\ < ^N\ 
II *(*i) |- VJV / 



1 1 
64 2vr' 



By symmetry, in this definition we can replace Vi(Xx) with 1*2(^1) • 
We use the notation W a for the vector valued process W a = {Wl, W 2 ), 
where W^ and W 2 are independent Wiener processes with marginal 
distribution W?{t) - W?{s) ~ jV(0, a 2 {t - s)) VO < s < t, Va = 1, 2. 



Theorem 3.1. Let Zm be the process defined in ffTTj) . Then for any 
< T < oo ; {Zjsr(t)}o<t<T converges to the two-dimensional Wiener 
process {W a (t)}o<t<T- Convergence is in distribution on the space of 
continuous functions C ([0, 7], M 2 ) equipped with the uniform topology. 

Then we will prove that {T^ 1 (t)} ie [o i T] converges in distribution to 
the function t. Combining these two results, we can show that Z^oT^ 1 
converges in distribution to W a . Observing that Z^oT^ 1 is the process 

1 rNt 



this implies our main theorem. 

Theorem 3.2. For any < T < 00, {^Ar(t)}o<r converges to the 
two-dimensional Wiener process {W ff (t)}o<t<t<T- Convergence is in 
distribution on the space of continuous functions C ([0, 7], M 2 ) equipped 
with the uniform topology. 

Finally, we will use the previous result to show that the rescaled 
solution of the Boltzmann equation converges to a diffusion. Let us 
denote with u N the two dimensional vector-valued measure defined as 



u N (t,k,x) :=u{Nt,k,{N In N) 1/2 x), Vt > 0, \/k e T 2 , Vi€ 



e. 



where u is solution of ([3]) in dimension two, with initial condition 

u(0,k,x) = u (k, (N\nN)- 1 / 2 x). 

Theorem 3.3. Assume that u E L 2 (R 2 x T 2 ;R 2 ) and fix T > 0. 
Then, Vt G (0,7], u N (t, •, •) converges weakly in 7 2 (R 2 x T 2 ;R 2 ) to 
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u(t, •), which solves the following diffusion equation 



(14) 



d t u(t,r) = -a 2 Au(t,r) 



u a (0, r ) = ~J2 I dk u o ( r > k ) Va e x > 2 ' Vr G 

2 ,3=1,2 ^ 



4. Sketch of the proof 
4.1. Theorem 13.11 Define the two-dimensional random vector 
(15) Vn := $(X„)-^(X n ), n G N„. 

We will denote with if)%, a = 1,2, the a-component of ^„. 



We decompose Zjv, defined in (TTTl) in two parts, i.e. Zjy = ^ + Zf 



where Va = 1, 2, Vt > 0, 



IJWJ-l 



Zft(t) = (iVlniV)- 1 / 2 £ e n ^l {enm>vw} 

n=0 



+ (iVlniV)- 1 /2 eL]v ^ rjvtj i {eLivij| ^ |> ^ } (iVt _ [Nt \) 

LJVtj-i 
Z£<(*) = (iVlniV)-^ £ e n Cl K |«|<^ } 

n=0 

+ (iVlniV)- 1 /2 eL]vtjV;r]vtj i {eLivtj| ^^ v/F} (iVt _ Li vtj) . 

p 

At first we will show that Z^ — ¥ when iV — > oo. It is enough to show 
that for every unitary vector A := (Ai, A2) 

XiZ 1 ^ + x 2 z 2 > 4 0, iv -»■ 00. 

This is stated in the next Lemma. 
Lemma 4.1. For every 5 > 



(16) 



lim P 

7V->-oo 



sup \X 1 Z 1 N > (t) + X 2 Z 2 N > (t)\ >5 

t£[0,T] 



VA G R 2 such that \X\ = 1. 

Proof. For every A G K 2 with |A| = 1 



P 



sup lAiZ^OO + AsZ^OOl >5 

«6[0,T] 



< P 



sup {|^ > (t)| + |^ > (t)|}>5 
te[o,T] 



0=1,2 



sup \Z%>{t)\> 5 - 
«e[o,T] z 
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For every t G [0, 7], Va = 1, 2 

|Zi5>(i)l - TfEf S e " ICI Vi«>^}' 

Then, by Chebyshev's inequality 



P 



2 
< - 



snp \Zf(t)\> 5 - 

*6[0,T] Z 



v n=0 

where in the last inequality we used the fact that \/n > 0, Va = 1,2 



E 



_ e «ICI 1 { e „|^|>^}J < Co^, 



as one can easily compute, using the upper bound for P m 021]) and the 
fact that \k\ 2 \ip a (k)\ is finite for every k G T 2 , Va = 1,2. 



as one 



U 



Let us consider Zf T . As first step, we will prove that for every unitary 



N ' 
vector A G M 2 , (Z£, A) := XiZjf + \ 2 Zj^ => W ff , where W a is a one 

dimensional Wiener process such that W a (t) — W a (s) ~ A/"(0, a 2 {t — s)). 

This is stated in the following proposition, the proof is postponed to 

the next section. 



Proposition 4.2. Fix T > 0. Then as N — >■ oo, /or even/ A G M 2 , 
wrf/i |A| = 1, (Z^,\) converges weakly to the one dimensional Wiener 
process W a . Convergence is in distribution on the space of continuous 
functions on [0, T] equipped with the uniform topology. 

Now we have to show that Z^ converges to W a . We follow the ap- 
proach of (27] (see the proof of Lemma 4) . The tightness of the sequence 
{Z^} N >x follows from the tightness of the sequence {{Z^-, A)}jv>i, for 
every unitary vector A. Thus we only have to prove the convergence of 
the finite dimensional distribution. In particular, we have to show the 
following: 

(i) Z<(t) - Z<{s) =► W ff (t) - W„(s), V < s < t < T; 
(ii) Zh(s) and (Z^(t) — Z^(s)) are independent, as N — > oo, 

v o < s < t < r. 

In order to verify the first condition, we observe that the convergence of 
the process (Z<(-), A) to W„{-) implies that ((Z<(s), A), (Z<(t), A)) => 
(W ff (s),W^(t)), for every s,t > 0. But (W a {s),W a {t)) has the same 
law of ((W a (s),X), (W^t), A}), then 

(Z<(t), A) - (Z<(s), A) => (W a (t), A) - (W a (s), A) 
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for all V < s < t < T, VA G M 2 with |A| = 1, and this implies (i). 

In order to verify condition (ii) it is sufficient to prove that Z^(s) and 
Ztf(t) — Z N (s) are asymptotically jointly Gaussian and uncorrelated. 
This is stated in the next Lemma. 

Lemma 4.3. For all X, fj, G M 2 

(17) 
(Z<(s), A) + ((Z<(t) - Z N (s)),») =► M (0, a 2 {\X\ 2 s + |/i| 2 (t -a)}), 

v o < s < t < r. 



We postpone the proof in section 15.21 

4.2. Proof of Theorem 13.21 Converge in probability of T^ 1 to the 
function x, where x(t) = t, in a compact [0, T], is proved as in [3], see 
Lemma 8.1 and Proposition 8.2. Then 

(Theorem 3.9 in Billingsley [6J) and therefore Z N o T^ 1 =>■ W a o x 
(Billingsley jB], Lemma pg. 151). 

4.3. Proof of Theorem 13.31 Given a complex vector valued function 

J on I 2 x T 2 , we define the Fourier transform on M? x T 2 

j(p,k)= f due- ip - u J{u,k). 
Jv? 

We choose a class of test-function J G S{R 2 x T 2 ; C 2 ) and we introduce 
the norm 

||J|U= / dp sup \J(p,k)\ 2 . 
Moreover, we denote with (J, F) the scalar product 

(J,F)= V / dp [ dkJ a (p,k)*F a (p,k), 

Q=1)2 JR2 JW 

which is well defined for every F G L 2 (M. 2 x T 2 ). We use a probabilis- 
tic representation of the solution of the rescaled Boltzmann equation, 
namely 

(J,u N (t)) 

= J2 dpdkJ a (p,kyE {a:k) [u (p,a {Nth K (Nt) )e- tp - YN{t) ], 

where E( a ,£;)[- ] is the expectation starting from the state (a, k), and 

We choose a sequence of real numbers {0n}n>i such that On — > oo 

for N | oo and ^ff= -)- 0. Since Vt G (0, T] 



VNlnN 



e -ip-Y N (t) _ e -ip.Y N (t-^t) 



< c -^^\p\r, 

VNlnN 
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with Cq finite, using Cauchy-Schwartz we get 



E 

a=l,2 



2 xT 2 



dp dkJaipi ky 



<cr 



Uo(p,a(m),K(Nt))(t 



i P -Y N (t) _ ~i P .Y N (t-^ty 



7 iV 



VN\nN 



|V U J|L \\uo 



e 



A \\ u 0\\L' 2 (R 2 xT 2 ) 



0. 



The measure n on {1, 2} x T 2 , given by Tr(a, dk) = ^dk, is invariant 
for the (reversible) process {(a(t),K(t)),t > 0} on ({1,2} x T 2 ). De- 
noting with U p (a t ,K t ) = u (p,a t ,K t ) — tt[u }(p), \/p e M 2 , Vt > 0, we 
have 



E 



(a,k) 



(u (p,a (Nt) ,K {Nt) ) - n[u }(p))e 



-i P -Y N (t- e -$-t) 



E 



(a,k) 



-iV-Y N {t-^t) g^ U p ( at _ eNt7 Kt _ dNt) 



where {S t }t>o is the semigroup associated to the generator Q.Thus, 
using Cauchy-Schwartz, 



(19) 



E 

0=1,2 



x E 



dp dkJ a (p, k)" 



! xT 2 



(a,k) 



< 2 \\J\ 



A 



(u (p, a (Nt) , K {Nt) ) - 7r[u ](p))( 



dp \\Se N t Mp\\l 2 -(t 2 )- 



-ip-Y N (t- H -$Lt) 



In order to prove that St converges exponentially in the L\ norm, we 
have to check that — (Cf, /)- > A ||/|| L % , with A strictly positive, for 

every two dimensional vector valued / G Lf (T 2 ) such that (/, 1)# = 
(spectral gap). This condition can be easily verified by observing that 

P(k,k') 

a:= sup ^zk 2 ^ 1 ^, 

fc,fc'eT 2 <P{k') 
thus 

- (Cf, /> # > (1 - o) / dk<f>(k) \f(k)\ 2 > 0, V/ e L\ (T 2 ) : / JL 1. 
This implies 
(20) / dp \\Sg Nt W P || L 2( T 2) < e ° ^ ||mo|Il2( R 2 xT 2) , 



with Aq > 0, then the r.h.s. of (TT9I) goes to zero for N — >■ oo, Vt > 0. 
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Finally, by Theorem G£2J E (a>fc) [exp{ip ■ Y N {t)}} -» exp{-(p?H-p|)o- 2 t/2}, 
Va = 1, 2, Vk G T 2 /{0}, then using (ITSjI . (fT9l. (120)1 we get 

(J,u^(t)> ^ V / o>^ .Up,kYn[u ](p)e-^ 2t / 2 , 

ViG (0,71. 

5. Details 

We start with some preliminary results on P m , the m— th convolution 
integral of P, the probability kernel defined in[HJ By direct computation 

(21) P m (k, dk') = 2 fl E sin2 ( 7r ^)4 r 3 rin 2 (^)dA/ 

E 7 =i sin (vrA; 7 ) ^ ^ 

where, Va, (3 G {1,2}, 

(22) 4^ = <W, 4T 1) = t am U Vm>l. 
Here a is a 2 x 2 real matrix with elements 

sin 4 (7r/ci) 



On =^22 = 2 / d£; 

JT 2 

O12 = Q2i — 2 dk 



T 2 "'" E a Sm 2 (7r£; a )' 

sin 2 (7rfci) sin 2 (7rA;2 



^T 2 Z)a Sil1 (^ 

Observe that the condition 



/ P m {k,dk') = 1 Vm> 1 



implies 



2 



(23) E A S = 1 ' V« = l,2, Vm>l, 

/3=1 

and thus 

(24) P m (A;, dife 7 ) < 2 ^ sin 2 (irk' )dk', \/k G T 2 , Vm > 1. 

,3=1,2 

5.1. Proof of Proposition [Q Fix A := (Ai,A 2 ) with A 2 + \\ = 
1. We will follow the strategy of Durrett and Resnick [8] to prove 
that (Ztf, A) := \iZlf~ + XiZ 2 ^ converges weakly to a Wiener process 
W c . They use a result of Freedman [13], pages 89-93, on martingale 
difference arrays with uniformly bounded variables. We start with the 
following 

Definition 5.1. A collection of random variables {£jv,i}, N > 1, i > 1 
and a - - fields J 7 ^, i>0,iV>lisa martingale difference array if 

(i) for all N > 1, J-N,i, i > is a nondecreasing sequence of cr-fields; 
(ii) for all N > 1, z > 1, £jv,j is -^jv.j measurable; 
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(iii) for all N > 1, E [Cn^I^na-i] = a.s. 
We introduce the following notations: 

(A '^' m):=Al ^^ 1 {-l^^} 

1 j e ^ 2 

ViV > 2, m > 0, and, for JV = 1, 

(A, *i >m ) = Aie m V' m l{e m [^i,[<i} + A 2 e m ^ m l{ em |^2 i |<i}, 

for all m > 0. 

For all iV > 1, m > 0, we denote with J r N,m the cr-field generated 
by {X , ..,X m } x {e , .., e m }, where {A" m } m > is the Markov chain with 
value in T 2 . Then {(A, ^ N,m) , J^n, m} N>i,m>i is a martingale difference 
array. In particular, condition (iii) of 15.11 can be easily checked using 
the explicit form of probability kernel P[k, dk'} . 

By definition, the variables (A, ^N,m) are uniformly bounded in m, 
i.e. for all AT > 1 | (A, ^N,m) | < £jv ? V«i > 0, where e^ = ,, 2 N if 
N > 2, and E\ = 2. In particular e^ I when A/" — >■ oo. 

For every N > 1, j > 1, let us define 

i 

(26) (A,^-> = ;jT(A,^v, m >, 

m=l 

J 

(27) (A,y^-) = ^[(A,^) 2 !^,^]. 

m=l 

We will prove in lemma l5~2l that P [lim^^A, Vnj) = oo] = 1, for all 
N > 1, i.e. the martingale difference array {(X,^ N,m) , J^n, m} N>i,m>o 
satisfies the hypotheses of Theorem 2.1 in [8]. Thus, setting 

JN,x(t) = su-p{j\(X,V N>j ) < t}, 

we get that (A, SjvjVa(-)) converges weakly as a sequence of random 
elements of D[0, 7] to a standard Wiener process W. 
Now let 4>N,\{t) = (A, Vjv.Livej)? V^ £ [0,7]. By definition 

In order to prove that 4>n,\ converges in probability to the function </> : 

4>{t) = cr 2 t, it suffices to show that (pN,\(t) — >■ a 2 t, Vt G [0,7], since <fi 
is continuous and <J)n,x is monotone. That will be proved in lemma I5T21 
Then 

(Billingsley [5], Theorem 3.9) and therefore 

(A, S NiJN J o(j) NiX =>• W o(f) 
(Billingsley [6], Lemma pg. 151). 
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Finally, 

(A,SV,ljv-j) = (A, SjvjatOM-))) ^ W<r> 

where convergence is in distribution on the space D[0, T] equipped with 
the Skorohod Ji-topology. 

The process (A, £#(£)) := 5^m=o (A, ^?N,m) converges also to W a . 
For every N > 2, (Z^,\) = XiZj^ + A2Z 2 / is the continuous func- 
tion defined by linear interpolation between its values (X,S^f(m/N)) 
at points m/N. The two sequences {(A, S^(t)), < £ < T} and 
{(■Zat^), A)0 < t < T} are asymptotically equivalent, i.e. if either 
converges in distribution as N — > oo, then so does the other. Conver- 
gence of {Zx, A) to W a is in distribution on the space of continuous 
functions equipped with the uniform topology. 

We conclude this subsection with the main Lemma. 

Lemma 5.2. For every N > 1, for every unitary vector AeK 2 , 



(28 



P 



lim (A, Vjvj) = oo 



1. 



Moreover, for every 5 > 0, for every unitary vector A 6 M 2 , 



(29) 



lim P 

JV-KX) 



(X,V NAm ) - a 2 9\ >5 



0. 



V0e [0,T]. 

Proof. Fix A G R 2 , with |A| 2 = 1. VJV > 2, we define /W : T 2 ->• M 2 



(30) 



/at (A;) = / cfe e~ 

'0 



x 



f ( v^ -2 ib a (k') 

U m dk>) { E Aa 7ivh7l 1{2 i^)i^ 



Using ([21]), we get /^(Ar) > C /iV, with < C < oo. Since 

/jv(X m ) = E [(#iv,m+i, A) 2 | 7-J , Vm > 

then, for all N > 1, (A, Vjvj) > j CoiV -1 which goes to infinity for 
j — > oo, a.s. 
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Now we focus on (1291) . By Chebychev inequality, for every N > 1 



P 



(\V N>lmi ) - aH\ >5 



< P 






n=l 

[Nt] 



-) 



>s- 



N 



(31) 



<^E E [( E [(^- A ) 2 i^-i]-^) 
^EE E (e[<^,„,a>v»-i]-^ 

=1 m^n L 

2 

(E[(*jv,m,A) 2 |J-m-l] "^ 



N n=l m^=n 
X 



where 5 N = 5 - N' 1 . By (JM]), we get 



(32) 



Co 



E [(<f> N , m ,\) 2 \Fm-l] = /v(X m _0 < -£, 



iV 



thus the first sum on the r.h.s. of (I3T1) is bounded by Sj^CxT/N, with 
C\ finite. Let us consider the second sum on the r.h.s. of (|3Tj) . For 
n > m 



E 



E [(* Ntn , A) 2 | F n -{\ E [(V N , m , A) 2 | Fm-i] 



E 



We set 



E [(^N^^lFm-i] E[E [<^jsr,«, A> 2 I^=k_i] |J" m -i] 



^- m (X m _i) := E E [(**,„, A> 2 |Jv.i] |J- m _i 



where, for every / > 1, iV > 1, the function g : T 2 — > R 2 is given by 



g l N (k)= / dk' P l (k,dk')f N (k'), 

Jt 2 



with /jy defined in®. By flMD and Q32J) we get 



(33) 



^(fc) < ^, Vk e T 2 , V/ > l. 
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We fix M, 1 < M < N and we get 



[Nt\ 

J2J2 E [ E [<**>»> A > 2 !^n-l] 1 [(*JV,m, A) 2 | J- m _J 
n=l m^=n 

M [Nt\ 

= 2^ ^ E[/ JV (x m _ 1 )^- m (x m _ 1 : 

m=l n=m+l 

|MJ m+M 

m=M+l n=m+l 
[Nt\ [Nt\ 

+ 2 J2 E E^A^)^™ (X m _!) 

m=M+l n=m+M+l 



L m— 1 , 



By fl33|) . the first and the second sum on the r.h.s. are bounded form 
above by CTM/N, with C finite. We denote by {iP m ~ x the convolution 
integral of the initial measure \x and the probability p m_1 . For every 
I > 1, 



E 



fN{Xm-l)9N\Xm-l) —^n fN{X m -i)g N (X m -i) 



+ / ^r-Vfc) - Tr(dfc)] f N (k)g l N (k) 



where the last term is bounded by CAT 2 J T2 \[iP m 1 (dk) — n(dk)\. 
Moreover, for every I > 1 



E„ 



fN{X m _x)g N {X, 



171— 1 , 



ir(dk)f N (k) / dk'P l (k,dk')f N (k') 
< ( / 7r(dk)f N (k)) +§ij T2 \^P m -\dk) - n(dk)\ 



We get 



^ ^E[E [(^ n ,A) 2 |7-n-l] E [<*JV,m,A) 2 |7- m _i] 

n=l m^n 

<L^*J(L^*J-i)(e. [(*jv,i,A> 2 ]) 2 
+ CT— + CT f \fxP M {dk) - 7r(dfc)| , 
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we C and C finite. In the same way one can prove that 

[Nt] 

J2^[(^N,n,X) 2 ] < [Nt\E n [(^,„,A) 2 ] 



71=1 



+ CT— + C'T I \fiP M (dk) - ir(dk)\ 



with some C, C finite, and finally we get 



P 



|(A,Viv iL jvt|) -v 2 t\>$ 



1 „~-M 
<— CT— 



+ -J- C'T I WP M (dk) - ii{dk)\ , 



where C, C are finite. (129|) is proved by sending M, N — y oo in such a 
way that M/N ->■ 0. 

D 

5.2. Proof of Lemma (14.3ft . We use the central limit theorem for 
martingale difference array ([9], Theorem 1; see also JTD], [H]) which 
states the follows: fix t > 0, and let {£,N,i, 3~N,i\N>i, i>o be a martingale 
difference array such that 

[Ntj 

(z) ^E[^ )i |7- JVl i-i]4ct, iVtoo; 

t=i 

Ljvtj 

(«) ^E[^l { |^ i | >e} |J- iV , i _ 1 ]4 0, JVtoo, Ve>0. 

t=i 

Then 

^6^=^(0,^). 

j=i 
By definition of Z<, VA G K 2 

(34) (A,Z<(i)> = (\,S N>m ) + (Nt- [Nt\)(\,* imi ), 

Vt G [0,7~], where (X,Sn,-) is defined in f[2T)j) . The rightmost term 
in (134")) goes to zero in probability by Chebyshev's inequality. We fix 
A, fj, G M 2 and < s < t < T, and we define the following array of 
variables: 

(\,* N ,i) if 0<2< [,Ns\ -1, 
(fi,V Nti ) if [Ns\ <i, VN>1. 



u ■ 



We denote with Tn,i the a- algebra generated by (X , ..., JQ) x (eo, ••, e 
V7V > 1, i > 0. Then {£jv,i, ^at^Iat^i, i>o is a martingale difference 
array. In particular, since \(u, ^N,i}\ < 2 (In N)~ l l 2 for every i > 1, for 
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every unitary vector v e M 2 , it follows that Vfe > 0, there exists N such 
that |£jv,i| < e i VAT > N, Vz > 1. Therefore condition (n) is satisfied. 

Moreover, with similar arguments of the proof of (T2TJ1) . one can prove 
that 

^ E [g^-F^-i] ^ v 2 \M 2 s + a 2 \fi\ 2 (t - s), 

i=i 

with a 2 defined in ffl3l). Thus 



i=l i=|_iVsJ i=l 

=>Af(0,a 2 {\\\ 2 s + \n\ 2 (t-s)}). 

6. An invariance principle for centered, bounded random 

variables 



In this section we present an alternative proof of Proposition 14.21 
We start with a CLT for arrays of centered, uniformly bounded random 
variables, based on the convergence of the moments to the moments of a 
normal distribution. Some asymptotic factorization conditions, holding 
on average, are required. Then we will use it to show that for every 
unitary vector A e M 2 , (\,Z<(t)) = XiZ l <{t) + X 2 Z 2 N < (t) => W ff {t), 

Vte[0T]. 

Proposition 6.1 (CLT). Let {X n ^ i = l,..,n,n > 1} be an array 
of centered random variables and suppose that exists e n 4- such that 
\X U} i\ < e n , for all n and i. Let S n = YH=i X n ,i- Then S n =>■ A/'(0, c), 
if the following conditions hold: 

(i) W > 1, for every sequence of positive integers {pi,--,pe} such 
that 3pj = 1, j e {!,.-,£} 

« — -, _ _ n\oa 

£ E[(X nili r...(X nM )^]^0 
(ii) W > 1 

n nfoo 

J2 E[(X nA ) 2 ...(X nM ) 2 ]--+c i 

Proof. The proof is based on the convergence of the moments of S n . 
Of course M[S n ] = 0, while for the second moment we have 

n n 

E [(Sn) 2 } =J2 E [(^n) 2 } +J2 E [Xi,nXj,n] ~> C, 
i=l tyj 

since the second sum goes to zero for condition (i). 
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Now let us compute the third moment: 

n n n 

E [(Sn) 3 } = J2 E [(^,™) 3 ]+3^E [(X hn ) 2 X J , n }+ Y, E [X hn X hn X Kn ] 

The last two sums go to zero for condition (i). For the first sum we 
have 



E E [( X ^) 3 ] 



i=i 



<J2 E KAn) 2 |An|] <^E E K**,") 2 ] ~ £ « C — >■ ° 
t=l i=l 

In the general case, the ra-th moment E [(£ n ) m ] is made up of terms 
of the form 

n 

A(pi,..,p t ) Yl E [(Xh,n) pi -(Xu,n) Pt ] , l<£<m 

with {pi, i — 1, ..,£} positive integers such that p\ + p 2 + .. + pi = m. 
Here A(pi, ..,pe) is the number of all possible partitions of m objects 
in £ subsets made up of p±, ..,pi objects. Since all sums containing a 
singleton (i.e. there is a p. t = 1) go asymptotically to zero, we consider 
just the cases with p, t > 2, Vz = 1,..,£. Observe that this implies in 
particular that £ < m/2. In this case 



J2 E[(x iun r...(x il!n y*] 



n 



<e™~ 2i E E [(X il , n ) 2 ...(X k , n ) 2 } 



r .m-2£J 



which goes to zero if £ ^ m/2. Therefore all odd moments are asymp- 
totically negligible, while for even moments asymptotically 

n 

E[(S n ) 2k ]~A k J2 ^[(X tl , n ) 2 ...(X tk , n ) 2 ]^A k c k , 

where A k is the number of all possible pairings of 2k objects, namely 

A k = (2k - l)(2k - 3) • • • 1 = (2k - 1)!! 

Finally 

F 17 <? "i m l ""^f J ^ m 

^ Y^n) j ■ -> y^ m _ ^jj cm/2 m euen5 

which are the moments of a Gaussian variable J\f(0, c). □ 

Let us consider the array of variables {(A, ^?N,m)iN > 2,ra > 0} 
defined in (j25p . (1151) . with A G K 2 unitary vector. We have 

LJvtJ— i 
(X,Z<(t))= J2 (\*N, m ) + (Nt-[Nt\(\,* N>m )), 

m=0 

Vi G [0,7~], ViV > 2, where the rightmost term goes to zero in proba- 
bility by Chebyshev's inequality. By definition, (A, ^N iTn ) < , 2 for 
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every m > 0, ViV > 2. Moreover, since ip(k) is an odd function, and 
the probability kernel P(k, dk') has a density which is even in both k 
and k', the array satisfies condition (i). In order to check condition 
(ii), we will use the following Lemma. 

Lemma 6.2. For every £ > 1, for every sequence (mi,....,m^) such 
that rrii > 0, rrii > 1, for every N > 2 

(35) E [(A,^ mi ) 2 ....(A,^v, mi+ .. +m ,) 2 ] < j± 

with Co finite, Vt G [0, T] . 
Proof. By definition 

E [(A, ^7V, mi ) 2 ....(A, ^N,m 1+ ..+m e ) 2 ] 

= / dz ie -*- / //P mi (*i)(A,*iv(A;i^i)) 2 / ... 

Jo Jt 2 J 

poo p 

x / dz m e~ Zm I P m *(/c m _i,/c m )(A,^7v(/cm,2: m )) 2 
Jo Jt 2 

<2 e (J°° dz e~ z J j(k)(\^ N (k,z))A , 

where in the last inequality we used (L241) . We conclude the proof by 
observing that 

Km TV I dze- z 7t(k)(\,y N (k,z)) 2 = a 2 , 

N^oo J Q J J2 



with a defined in (JT3]). D 

We observe that 

J2 E[(\,* N>il ) 2 ....(\,* N)it ) 2 ] 



e{o,..,LAftj-i} 



i] E E E [< A > *iV, mi ) 2 -.(A, ^, mi+ .. +m ,) 2 ] • 



mi>0 m2,..,m^>l 

mi+...+m£<LAftJ— 1 



We split the sum on m\ in two part, namely ^2 mi=0 + ^2 m >mi w hh 
< M < [Nt\ - 1. Using ([25]) and the relation 

lim V N~ k = 1 

AT->oo ^—' fc! 

171 ^ ,„„,mi, >1 

we get that for every £ > 1, TV > 2, V£ G [0, T] 

A/-1 

^ ! E E E [< A > ^, mi ) 2 ....(A, ^, mi+ .. +m ,) 2 ] 



mi=0 m 2 ,..,m^>l 

m 1 + ...+m i <\_Nt\—l 



- u n 
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By repeating this procedure for all the sums, we have 

(36) 

J2 E[(\,* NM ) 2 ....(\,* Ntil ) 2 ] 



6{0,..,LiVtj-l} 



£\ 



J2 E[(A,^ imi ) 2 ....(A,^ mi+ .. +m ,) 2 ] +S £ (M,N), 



mi,.. ,mg >_M 
^ 1 + ... + m e <lNt]~l 



with S £ (M, N) < C t T t - 1 M/N, W > 1. 
Observe that for every m > 2 



P m (k,dk')(\^ N (k',z)) 2 = / 7i(dk')(X,^ N (k',z)y 

T 2 Jj 2 



+ 



T 2 



P m - 1 {k,dk)-ir{dk) / P(k } dk'){X,^ N (k' \z))\ 

T 2 



where, using (124ft. 



sup 

fceT 2 </t 2 



P m -\k,dk)-7r{dk) I P(k,dk')(X,^ N (k',z)Y 

T 2 



.Co 
- ~W sup , 

iV fceT 2 J J 2 



P m - 1 (k,dk)-7r(dk) 



Thus, thanks to fl35l) . for every (mi, .., m^) with m^ > M, i — 1, . 

E [(A, ^jv, mi ) 2 ....(A, ^jv, mi +..+ m ,) 2 ] 

(37) / /-oo /• _ \< 

= ( / tfee" 2 / Tr(dk')(\,^ N (k',z)) 2 \ + e t (M,N), 

where 

e^(M,iV)<£-^ sup sup/ P m (k,d~k)-n(dk) 

N m>M-l fceT 2 Jt 2 

Finally, by ( 136]) and ( 137]) we get 

]T E[(A,^ 1 ) 2 ....(A,^, J ,) 2 ] 



e{o,..,LJV*J-i} 



£\ 



J2 (E^[(A,^,!) 2 ]) £ + ^(M,iV), 



mi ,.. ,mf>Al 
mi + ...+mg<[Nt\— 1 



where 
(38) 



M 



K e (M,N) < CfT L [ — + sup sup 

-/V m >M-lfceT 2 JT 2 



P m (k,dk)-7l(dk) 
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In the limit M, N ->• oo such that % ^ 0, TZ e (M, N) ->■ and 



t\ J2 (E 7r [(A,f JV>1 ) 2 ])'->(cr 2 )V 



.2^ 

zt + .^ + mg^lNtl-l 



with a defined in (|T3|) . Thus the array of variables {(A, \I/jv, m ),iV > 
2, m > 0} satisfies also condition (iz), and we get 

[Nt\-1 

S N {t) := Y. ^ *"•»> ^ ■ A/r (°> ^ *)' 

n=0 

Vt G [0, 7^, VA G M 2 such that |A| = 1. 

We can easily adapt the proof and show that V0 < s < t < T 

S N {t) - S N {s) -»■ jV(0, a 2 (t - s)). 

In order to prove the convergence of the finite dimensional marginal 
to the Wiener process W a , we have to show that Vn > 2, for ev- 
ery partition < ti < ... < t n < T the variables S^iti), S^fe) — 
Sjsr(ti),..,SN(t n ) — Sjv(^n-i) are asymptotically jointly Gaussian and 
uncorrelated. This is stated in the next Lemma. 

Lemma 6.3. For every n > 1, V a(n) := (a%, ..,a n ) G lR n such that 
\<x(n)\ = 1 

n / n 

(39) J2 <Xk{S N (tk) ~ 5jv(i fe -i)) ^Afio, a 2 J2 <4(h - **- 

fe=i V fe=i 

V0 = t < ti < .. < t n < T. 

Proof. The case n — 1 is proved. Let us consider the case n — 2. Fixed 
(ai,a 2 ) G M 2 , with a 2 + a 2 = 1, we consider the following array of 
variables 

t,N,m = {aillmKlm-LJ-l} + «2l{m>LiVtiJ}) (\ ®N,m), VAT > 2, Vm > 0, 

which are uniformly bounded by -t= and satisfy condition (i). Let us 
define, Vt > 0, m > 0, N > 2, 

0-N,m(t) := C*ll{m<LiVtJ-l} + a 2 l{ m >LiVtJ} > 

which is uniformly bounded by 1. In order to check condition (ii), we 
repeat the steps done for S>N(t) and we get 



E E [&*■■■&*] 



e{o,..,LJV*2j-i} 



= *! E ajv,u(*i) 2 -ajv,v(*i) 2 (^ [(A,^,i> 2 ])' 

0<i 1 <..<^<LiVt2J- 1 

+ K t (M,N), 
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with 1Zi(M,N) the same of (|3"8|) . By direct computation 



'J 

0<i 1 <..<i f <LiVt2j-l 



E^ ! E («i) 2fc E w 



fc=0 l<i!<..<i fc <LiVtlJ lNt 1 j<i k+1 <..<i e <lNt 2 i 

then using 



E N ~ k N ^° lv Ni ( E ^ [< A '^> 2 ])' ^ (* 2 )'> 



l<i 1 <..<i fc <iV 



with a defined in ffT3|) . we get that condition (ii) is satisfied, i.e. 

&, E E [^-4J 



iV-*oo 

e{o,..,Livt 2 j-i} 



^ k\{l - k)\ al tla2 {h ~ tl] 



^Yialh + aUh-h)]' 



thus 



L-MfcaJ-l 
OLlSfffa) + a 2 [S N (t2) ~ S N (t x )] = ^ ^N,m 

->tf(0,(o 2 )[alt 1 + <4(t 2 -t 1 )]). 

The proof can be repeated for n > 3, in that case we find the multino- 
mial formula for a polynomial with n terms to the power £. 

□ 
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